Coulomb-Enhanced Spin-Orbit Splitting: The Missing Piece in the Sr 2 Rh04 Puzzle 
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The outstanding discrepancy between the measured and calculated (local-density approximation) 
Fermi surfaces in the well-characterized, paramagnetic Fermi liquid S^RhCU is resolved by including 
the spin-orbit coupling and Coulomb repulsion. This results in an effective spin-orbit coupling 
constant enhanced 2.15 times over the bare value. A simple formalism allows discussion of other 
systems. For Sr2Rh04, the experimental specific-heat and mass enhancements are found to be 2.2. 

PACS numbers: 71.18.+y, 71.20.-b, 71.30.+h 



Since the discoveries of high-temperature supercon- 
ductivity and colossal magnetoresistance in Mott insula- 
tors made metallic by hole-doping, transition-metal ox- 
ides have remained at the forefront of research. Their 
many lattice and electronic (orbital, charge, and spin) 
degrees of freedom are coupled by effective interactions 
(electron-phonon, hopping, t, Coulomb repulsion, U, and 
Hunds-rule coupling, J), and when some of these are of 
similar magnitude, competing phases may exist in the 
region of controllable compositions, fields, and tempera- 
tures. The interactions tend to remove low-energy de- 
grees of freedom, e.g. to reduce the metallicity. This 
rarely happens by merely shifting spectral weight from 
a quasiparticle band into incoherent Hubbard bands, 
as in the U/i-driven metal-insulator transition for the 
single-band Hubbard model, but is usually assisted by 
lattice distortions which break the degeneracy of low- 
energy orbitals and split the corresponding quasiparti- 
cle -or partly incoherent- bands away from the chem- 
ical potential. According to recent calculations using 
the local density-functional plus dynamical mean-field 
approximation (LDA+DMFT), such Coulomb-enhanced 
crystal-field splitting seems to be the mechanism trigger- 
ing the expansion-induced metal-insulator transition in 
undoped LaMnC"3 [l| and in V2O3 (2j, long considered 
the prototype Mott transition. The low-temperature, 
antifcrromagnetically-ordered, insulating phase of V2O3 
is well described [3j in the LDA+U static mean-field 
approximation, which yields the configuration t 2 2 —> 
e o a ig- Although this approximation exaggerates the 
tendency towards symmetry breaking, it does give a rea- 
sonable description of the shape of the Fermi surface (FS) 
on the metallic side of the transition [l|, [2| . 

When going from 3d to Ad transition-metal oxides, the 
larger extent of the Ad orbitals cause the hopping, t , and 
the coupling to the lattice to increase, and U and J to 
decrease. This is reflected in the rich electronic proper- 
ties of e.g. the t 2 4 g ruthenates in the Ruddlesden-Popper 
series (Cai-xSr^+iRu^C^+i 0, H, @, 0, @| ■ Here, the 
end- members (y—1 and i/=co) have the same structures 
as respectively La2Cu04 (2D K^NiF^type) and LaMn03 
(3D perovskite). The relatively small size and strong 
covalency of the Ca ions cause the Ru06 octahedra to 



rotate and tilt. The resulting misalignment of the Ru 
ti g Wannier orbitals (WOs) reduces the hopping between 
them, and so does the deformation of the WOs caused 
by Ca-0-i2g covalency [3, [10J. As a result, the 2D ma- 
terials with x < 0.1 are Mott-insulators. Ca2Ru04 is 
insulating below 360 K, has orbital order with flat octa- 
hedra below 260 K [ll(, and is antiferromagnetic below 
110 K with configuration xy^ xz \ y Z \ according to the 
LDA+U 12j,|l3]. Moderate pressure induces a first order- 
transition to a metallic state with reduced tilt and rota- 
tion, and with ferromagnetic order at low temperature 
l4 llol 1 161 ] . For x > 0.1 the materials are metallic and 
exhibit metamagnetism coupled to structural distortions 
as long as x < 0.25 17|,ll8(. The properties of the ruthen- 
ates with x > 0.1 seem to be well described [13, LL3, |20j, 
SI, [H, ES Q in the spin-unrestricted LDA (LSDA), a 



parameter-free approximation which essentially neglects 
U, and substitutes J by the Stoner exchange coupling. 
For Cai.5Sr .5RuO4 at 40 K, angle-resolved photoemis- 
sion spectroscopy (ARPES) gives a FS which -after pay- 
ing due attention to surface reconstruction- is found 25 1 
to have neither orbital nor spin-polarization, and to be 
in good agreement with the LDA. Finally, stoichiomet- 
ric Sr2Ru04 is tetragonal and becomes superconducting 
below 1 K, presumably with spin-triplet p-waye pairing 
[H H3- Both dHvA (ill and ARPES 



o-waye pe 
H E| mea- 
surements show that, at low temperature, Si~2Ru04 is 
a nearly 2D Fermi liquid whose FS agrees well with the 
LDA 13, [20j and a mass-enhancement of about 3. 

In view of this decreasing strength of the Coulomb 
correlations, it therefore came as a surprise when dHvA 
and ARPES (HSHIl] at ~10K showed Sr 2 Rh0 4 , also 
a paramagnetic Fermi liquid with similar structure and 
electron-electron interactions but one more electron, to 
have a FS in substantial disagreement with the LDA. 
The initial surprise was that the experimental FS has no 
xy sheet, but it was soon realized 
K2MF4 structure the 2D xy and x" — y 
broad that theyoverlap at the d 5 Fermi level, and that 
the observed [35j relaxation by alternating rotations of 
the neighboring, corner-sharing RhOg octahedra around 
their z-axes, is such as to gap those two bands. As a 
consequence, only the equivalent ID xz and yz bands, 



34j that in the 
2 bands are so 



TB+'SO' 
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FIG. 1: FS of Sr 2 Rh0 4 in the fc z =0 plane of the large Bril- 
louin zone centered at (k x ,k y ) — (0,0) = T, with corners at 
(tv, n) = M, and edge-midpoints at (n, 0) = X, in units of 
the inverse Rh-Rh nearest-neighbor distance. Grey: ARPES 
[3ll. [331] . the same FS in all six pictures, but unfolded in the 
top panels. Black: Theory using six different approximations. 
The top panels results from the analytical expressions and pa- 
rameter values given in the text and employs one xz and one 
yz orbital per cell. The remaining pictures result from all- 
orbital LDAJ36J] calculations for the proper crystal structure 
with 4 formula units per cell [351 ]. The conventional Fermi- 
sheet notation is given in the (LDA+U+SO) picture. 



which hardly hybridize with the other d bands nor with 
each other, remain at the Fermi level with the single t2 g 
hole distributed equally between them. However, also 
the LDA FS calculated for the proper structure deviates 
substantially from the experimental FS [3l|, [32], |33j . This 
discrepancy clearly seen in Fig.Q](LDA) is disturbing be- 
cause there is no experimental indication of any further 
distortion. Hence, Coulomb-enhanced crystal-field split- 
ting can not be the solution to this puzzle. 

In this Letter we shall argue that the clue is Coulomb- 
enhanced spin- orbit (SO) splitting. We begin by demon- 




FIG. 2: Theoretical band structure of S^RhCU in the six 
different approximations. See Fig. 1. 



strating that the experimental FS can be perfectly fitted 
and the low-energy band structure explained using sim- 
ple, analytical theory. These results are illustrated in 
the top panels of Fig.s Q] and O We then perform ab 
initio LDA calculations [36] including the SO coupling 
(LDA+SO), and later also on-site Coulomb effects in the 
LDA+U approximation. The detailed results in the 2nd 
and 3rd panels of Fig.s Q] and [2] prove the soundness of 
our simple theory. 

In the Hilbert space of only the xz and yz WOs, the 
2D translational symmetry is a quadratic lattice with one 
Rh per cell; the relative rotation of the nearest-neighbor 
RhOg octahedra merely modifies the hopping integrals 
[9(. In the z-direction the hopping is small and not de- 
tected in ARPES. Neglecting it, as well as the small 
rotation-induced hopping between the xz and yz bands, 
the electron dispersion is [3, |38[ : 



XZ — ^^7T COS fc'j- 



2t$ cos k v 



(i) 



and equivalently for Ey Z . Here, Ef is the position of the 
Fermi level with respect to the energy at f M/2 and signs 
have been chosen such that the hopping integrals are pos- 
itive: t v ~ fi dlT / (ep — e p ) because this hop is mainly 
via the Op z orbital, and < ts ^C t n because this hop 
is direct dd6. These tight-binding bands are shown in 
Fig.[2l(TB). In the space of the (xz 1, xz J., yz T, yz I) 



WOs, the eigenfunctions of the SO coupling are: 

X m] =\ = ( xz + iV z ) T and X-§ = {xz - iyz) \ 
with eigenvalue \C, (> 0) , and 

XI = (xz + iyz) I and x_i = (xz — iyz) f 

with eigenvalue — fC SO coupling thus splits the degen- 
eracy of the xz and yz bands along TM by £ as seen in 
Fig.[2](TB+'SO'). Since the structure has inversion sym- 
metry, all bands remain doubly degenerate. The band 
structure is given by: 






±Vfe-^) +c s 



(2) 



and the FS, £ k =0, is determined by merely 3 param- 
eters, t^/ep, ts/eF, and C/ £ f, which may be obtained 
by fitting to 3 FS dimensions, e.g. the two intersections 
along f M and the one along f X. Fitting in this way to 
the ARPES FS [331 ] . perfect agreement with the entire 
FS is obtained as seen in Fig^TB+'SO')- The area 
of the f-centered electron sheet minus that of the M- 
centered hole sheet equals half the Brillouin-zone area. 
As seen in Fig.[T](LDA+SO), the agreement with ARPES 
is less good for the ab initio relativistic LDA FS: the 
SO splitting along TM is too small. In fact, fitting the 
LDA+SO FS to our analytical expression, yields essen- 
tially the same values for t^jep and ts/ep, but Q/ef is 
smaller than Carpes/^f by the large factor 2.15. The 
LDA+SO calculation (Fig.gJ) yields: C = 0.13 eV, a 
value which is smaller than the 0.16eV obtained for el- 
emental fee Rh [39] due to the O p z tails of the rhodate 
WOs. Letting £ set the energy scale of bare bands, the 
TB+'SO' parameters reproducing the ARPES FS are: 
Carpes = 2.15 x 0.13 eV = 0.28 eV, e F = 0.260 eV, 
U = 0.185 eV and t s = 0.039 eV. The TB results in 
Fig.sQ] and [Hare obtained with ( = 0. 

In order to explain the origin of the SO -enhancement, 
we need to add the on-site Coulomb repulsion, H", and 
thus consider the 2-band Hubbard Hamiltonian: 



H = 



k fi—xz 7 yz a 



e k c k f c k 

c /i°/icr "-Tier 



E 



H so 



m 



(3) 



e k given by Eq. |T]) and n = (n x ,n y ) runs over the 
quadratic lattice. For simplicity, we drop the superscripts 
n in the following. The SO interaction is: 



H. 



so 



C 



§(' 






(4) 



where n m . is an electron-number operator and p is the 

|=i, i.e.: the SO po- 



diffcrcncc between 



if — tj and 



larization. The Coulomb repulsion is U between two elec- 
trons with the same orbital, U-2J between electrons with 



different orbitals and different spins, and U-3J between 
electrons with different orbitals: 



Hr 



U ?i3rii + n_3n_i 

\ 2 2 2 2 

+ (U — 2 J) ( fi3h_ 3 + Tii ti_i 

v ' \ 2 2 2 2 

+ (U — 3 J) ( 713 71_I + 7l_3 7ll 
v ' \ 2 2 2 2 



Spin-spin correlations have been neglected here and there 
is no double-counting correction because we consider 
equivalent orbitals. Since Sr 2 Rh04 is paramagnetic at 



low temperature, we can set 713 = ti_ 



(?i — p) and 



7ii = h_i = j (h + p) with 7i being the number of elec- 
trons in the two doubly-degenerate bands, and get: 



3U-5J 2 
He, ~ nT 



U-J 



U-J 



-pp. 



The last expression is the mean-field approximation, 
which grouped together with Hso © yields a one- 
electron Hamiltonian with SO-coupled bands ([2]), but 
with £ substituted by £ e // = (+ \ (U — J)p, and p deter- 
mined self-consistently. The polarization function, p (£) , 
may by found from the polarization of each Bloch state, 



'±;l 



dC/2 



±c 



y+e 



as follows from lst-order perturbation theory. As a re- 
sult: 



P(0 = 2C 



o Jo 



CLfcx CLr^ii 



(e* xz -e*y + C 2 



where the factor 2 takes the double degeneracy of each 
band into account and we have used that 9 (— £ k ) — 
6>(-e + ) = 6»(e + )_6»(-£ k ) . The integral is over the 
area between the f-centered electron sheet and the M- 
centered hole sheet. For £ large, p (£) — ► 4 — n, where ti 
is the number of electrons in the xz and yz bands, and 
for Q small, p (£) — > xoC with the susceptibility 



Xo 



x dk x dk y 9 (e k ) 9 (-e k ) 



1 



'XZ 



2 (U - U) ' 

(5) 
which turns out to be independent of n. The last ap- 
proximation results from integrating over the rectangu- 
lar £,5=0 area. The self-consistency condition, C, e ff — 
C + h(U — J)p (Ce//) , thus yields a Coulomb-enhanced 
SO coupling, which for £ e // in the linear range of the 
polarization function is given by: 



Ce// 
C 



1 



U-J 



-Xo 



1 



U-J 



4 (t„ - 1 5 ) 



(G) 



Inserting Ce///C = 2.15 and t^ — ts = 0.146 eV, we 
get: U — J = 0.3ey and p(( e ff) = 0.97, whereas us- 
ing the proper polarization function, which saturates at 



p(oo) =1, yields: U — J = 0.5eV\ This is a reasonable 
value for a Ad WO spreading onto the oxygen sites. 

In order to substantiate this simple picture, we perform 
all-orbital relativistic LDA+U calculations [36J with U- 
J adjusted such as to give the best agreement with the 
ARPES FS. Since U and J in such calculations do not re- 
fer to proper orbitals, but to d-waves truncated and nor- 
malized inside atomic (LMTO) or muffin-tin (LAPW) 
spheres, the values of the parameters depend on the 
sphere size and are generally larger than for the more dif- 
fuse WOs [11. As is obvious from Fig.[j](LDA+U+SO), 
this agreement is even more perfect than for TB+'SO'; 
now even the small observed gaps [37( induced by the 
AF rotations are reproduced. The LDA+U+SO bands 
in Fig. [3 are well reproduced in the range from 0.15 eV 
below- to 0.5 eV above the Fermi level by the TB+'SO' 
bands folded into the BZ/2 with corners at X. Features 
not reproduced are the tiny splittings due to in-plane xz- 
xy hopping and out-of-plane hoppings neglected in our 
TB model. The agreement between the TB and the LDA 
calculation is less satisfactory, first of all because with- 
out SO-quenching the in-plane xz-xy hopping produces 
a splitting along TM, and secondly because the rotation- 
induced xy-{x 2 — y 2 ) gapping is not complete without 
SO coupling. In fact, it takes the LDA+U+SO to push 
the lower edges of the xy- [x 2 — y 2 J gap to — 0.16eV 
along f M, and even deeper along f X, locations close to 
those observed with ARPES [32J. Note finally, that the 
LDA+U alone, without SO coupling, brings little im- 
provement compared with the LDA. 

Having obtained perfect agreement with the ARPES 
FS, we use the LDA+U+SO to calculate cyclotron 
masses and compare with those obtained from the dHvA 
measurements [32J. The resulting enhancements are 
toarpes/ttjlda+u+so = 2.1 (a) , 2.1 (/3 M ) , and 2.3 (/3 X ) 
Consistently herewith the density of states at the Fermi 
level yields the electronic specific-heat enhancement 
7cxp/7LDA+u+so = 2.2. These many-body enhancements 
are smaller than those (~ 3) in Si^RuOj. 

We have finally performed relativistic LDA+U calcula- 
tions also for S^RuOj. and find the agreement with the 
experimental FS 28l.l29l.l30l] to improve from good to per- 
fect. Inclusion of SO+U in the LDA reduces the area of 
the a-pocket from 13 to 11% of the BZ-area in Sr 2 Ru04, 
as compared with a 24 to 6% reduction in Sr 2 Rh04. The 
stronger SO effects in the latter oxide are caused by £ 
being 20% larger [39j, due to the larger mass, and by the 
Coulomb enhancement ((6]) being larger due to the reduc- 
tion of tn caused by rotation. Other factors are similar 
in the two oxides because the extra hole in Sr2Ru04 is 
accomodated in the xy band, which hardly couples with 
the xz and yz bands. It is conceivable that reduction of 
tit caused by Ca-induced rotations could make Coulomb- 
enhanced SO coupling of the xz and yz bands important 
in other ruthenates and rhodates not driven magnetic by 
the i-reduction. 



In conclusion, resolution of the Sr2Rh04 puzzle has 
tought us that although usually neglected in 4<i-oxides, 
the spin-orbit coupling belongs to the list of competing 
interactions which cause the rich physics of these mate- 
rials [|(J. 
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